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] .  INTRODUCTION 


In  the  wave  propagation  in  randomly  fluctuating  ocean,  most  researchers  arc 
concerned  with  the  time- independent  random  media  [1],  i.e.  the  local  sound 
speed  varies  independently  of  time.  For  short  range  propagation,  this  is 
perhaps  not  a  serious  omission.  However,  in  the  long-range  transmission,  the 
temporal  fluctuation  of  media  becomes  important  and  should  be  taken  into 
account . 

The  main  objective  of  our  work  is  to  investigate  the  statistical  laws  in 

the  acoustic  wave  propagation  through  random  media,  such  as  the  ocean.  We 

introduce  the  method  of  progressing  waves  for  the  time -dependent  random  media. 

If  the  fluctuation  is  weak,  this  is  the  case  in  ocean,  we  are  interested  in 

the  accumulative  effect  over  a  long  time.  This  consideration  leads  to  the 

2 

diffusions  limit:  c  -*■  0,  t  -*•  °°  with  t  =  z  t  fixed,  where  e  is  the  scale  of 
random  fluctuation.  Our  major  result  is  that,  in  the  diffusion  limit,  the 
statistics  of  the  phase  and  the  amplitude  for  each  wave  component  may  be 
computed  easily  according  to  the  probability  law  of  a  simple  diffusion 
process  with  the  variance  parameter  o  depending  on  the  covariance  function, 
the  invariance  principle  of  Donsker  in  the  classical  situation. 

In  Section  2,  we  develop  the  basic  theory  via  a  simple  model  equation. 

The  first-order  random  PDE  (partial  differential  equation)  is  to  be  solved  by 
the  method  of  characteristics.  Then  we  apply  the  Khasiminski 's  [2,3]  limit 
theorem  to  random  ordinary  DE  for  the  characteristics.  Thereby,  in  the 
diffusion  limit,  we  can  determine  the  fluctuations  of  phase  and  amplitude  by 
the  limiting  diffusion  process.  In  the  next  section,  the  system  governs  the 
acoustic  waves  in  one  dimension  is  treated.  There  we  introduce  the  method  of 
progressing  wave  solution.  This  together  with  the  method  of  characteristics 
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enable  us  to  resolve  the  question  of  statistical  fluctuations  of  waves, 
finally,  in  Section  4,  the  results  in  one  dimension  will  be  extended  to  the 
physically  important,  three-dimensional  problem. 

2.  A  SIMPLE  MODEL  EQUATION 

To  illustrate  the  basic  ideas,  we  consider  the  following  simple,  first-order 
partial  differential  equation  with  random  coefficient 

(2.1)  ♦  v(t,x)  ^L]  u  =  ne(t,x,u))u 

(2.2)  U|t=0  =  a(x) , 

where  a,  V  are  the  given  functions  and  r)£  is  a  random  function  depending  on  a 
small  parameter  c  >  0. 

By  the  method  of  characteristics  in  partial  differential  equations,  we 
introduce  the  characteristic  equation 

(2.3)  ^  =  V(x,y)  ,  Y(t)  =  x  , 
whose  solution  is  a  characteristic  curve  T  through  (t,x) 

(2.4)  y  =  C(s;t,x)  ,  0  <  s  £  t  . 

Along  T  the  system  (2.1)  and  (2.2)  becomes 

(2.5)  d*  u[s ,£(s;t ,X)]  =  n£[s,C(s;t,x,w)]U[s,?(s;t,x)]  . 

(2.6)  U[0,C(0;t,x)]  =  a[5(0;t,x)]  . 

The  solution  of  the  above  simple  ordinary  differential  equation  is 

(2.7)  U(s;t ,x,oj)  =  a[£(t ,x)]exp  {/*ne[a;£(o;t,x)]do}  , 


where  C(t,x)  =  5(0;t,x)  . 


From  (2.7)  the  solution  of  (2.1)-(2.2)  is  given  by 

(2.8)  U(t,x,w)  =  U(0;t,x,u))  =  a[£  (t  ,x)  ]exp{0e  (t  ,x,u)) }  , 

where  the  amplitude  a[^(t,x)]  is  deterministic,  while  the  phase  function 
(t  ,x,(i>)  =  Oe(s;t  ,x,u>)  |  s_q  with 

(2.9)  ee(s,t,x,u ))  =  /^ne[o,5(o;t,x) ,w]do  , 
is  random. 

For  fixed  (t,x),  the  dependence  on  which  will  be  suppressed,  (2.9)  takes 
the  form 

(2.10)  ej(<u)  =  /^ne(s,ai)ds  . 

In  studying  the  statistical  law  of  the  solution  (2.8),  one  naturally  raises  the 
question,  as  t  1  0,  (t-s)  t  °°,  whether  a  family  of  processes  {0E,  t  >.  0}  con¬ 
verges  weakly  to  a  certain  limi.  distrib  on  0t(w),  under  appropriate  scaling. 
The  weak  convergence  of  to  6  will  be  denoted  by 

(2.11)  .  0E  =>  0  . 

The  expression  (2.10)  is  in  the  form  of  the  "sum"  of  random  variables.  If 
the  necessary  moments  exists  and  the  process  ne(t,u>)  is  asymptotically  inde¬ 
pendent,  then  the  normalized  {0E  -  E  0E}/a{0E}  =  0E  *=>  0^ 6  V  (0,1),  the  normal 

distribution  with  mean  zero,  variance  one.  It  follows  the  phase  fluctuation 

2 

is  asymptotically  normal  with  mean  E0t  and  variance  o  {0^.}.  Thus  the  wave 
function  satisfies  the  log-normal  law. 

Another  interesting  case  is  the  diffusion  limit.  This  limit  is  important 
because  it  renders  the  problem  computable  and  permits  generalization  to 
higher-dimensional  situations.  Suppose  is  such  that 

(2.12)  0e  »  n0  +  en 

where  n0  is  deterministic,  Ef\  =  0  and 


(2.13) 


E(n(t)n(s)}  ■  a(t,s) . 
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Assume  the  limit 

t  +T  t  +T 

(2.14)  o  =  lim  ■=  J  f  a(t,s)dtds 

T-*°  1  co  co 

exists,  independent  of  t^,  and  the  process  is  asymptotically  independent  (to 

£  £  2 

be  made  precise  later).  Then  the  rescaled  phase  function  0  (t)  =  0  (t/t:  ), 
for  fixed  t  =  ezt,  we  have 

0C(r)  ->  G°(t)  , 

which  is  a  Wiener  process  with  the  variance  parameter  o.  This  is  a  consequence 
of  the  well-known  Khasiminski  theorem  [3]. 

Now  we  consider  a  variant  of  (2.1): 

(2. 15)  [~  +  ne(t,x,w)  -~]\}  =  b(t,x)  u, 

(2.16)  U|t=Q  =  a(x)  . 

In  contrast  with  the  previous  equation  (2.1),  the  characteristic  equation 
,  c 

(2.17)  =  nG(s,y\w)  ,  yc(t)  =  x  , 

is  stochastic.  The  limit  theorem  for  the  stochastic  differential  equation 

(2.17)  was  treated  by  Khasiminski  as  mentioned  before, 
bet 

(2.18)  ne(t,x,u))  =  nQ  +  en(t,x,uj) 

where  is  a  constant  (in  general,  may  be  a  function  x,t).  The  random  field 
satisfies 

(i) 

(2.19)  En  =  0 

(ii)  n(t,x,u))  is  (a.e.)  continuous,  and  its  partial  derivatives  up  to 
second-order  are  uniformly  bounded  in  [0,T]  *  F  . 

(ii  i) 

T 

(2.20)  o2  =  li*  i  J/En(t0«-T-t,  x_-ri0t)n(t0+T-s,x0-n0s)dt 

T-*®  0 
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which  exists  independently  of  t^  >  0  and  x^. 

(iv)  The  a-field  F*  is  increasing  in  [t,s]  such  that  n(t,x,u>)  is  measurable 
with  respect  to  FE,  and  for  all  A  in  BE  ^t+T’  we  ^ave 

(2.21)  | P (A  B)-  P(A)P(B) |  <  6(T)P(A)  . 

where  T^3(T)  I  0  as  T  t  ®  . 

Then  we  can  apply  Theorem  K  (Khasiminski ' s  theorem)  to  get 


2  2 

Lemma  2.1:  As  £  t  0,  (t-s)  +  °°,  with  t  =  e  t  and  o  =  e  s  fixed,  the  solution 
(2.17)  converges  weakly  to  a  Markov  diffusion  about  the  mean  characteristic, 
that  is 

(2.22)  YE(t)  =>  f.0(x,t-s)  ♦  C(t-a)  =  Yg(t) 

where  F0(x,t)  =  (x-nQt),  £(0,oj)  =  0, 


and  £(t,oj)  is  a  Markov  diffusion  process  with  mean  zero  and  the  diffusion  coeffi- 
2 

cient  o  defined  by  (2.20)#. 

Remark:  In  fact  the  above  lemma  says  that  the  limit  process  will  satisfy  the 
I to  equation 

(2.23)  dys  =  Hq  ds  +  eadw(t-s)  ,  y  =  x. 

Let  Tt(u))  be  the  random  characteristic  curve 

(2.24)  ye=  YE(t,ui),  Oisit  . 

Similar  to  (2.8),  along  rt(u)),  the  system  (2.1S)-(2.16)  may  be  integrated  to 
give 

(2.25)  U£ (t,x)  =  a[Y^(t)]exp{<Kt,yE)}, 

where 

(2.26)  4>(t,ye)  -  /5b[s,YE(t))ds  . 
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To  ensure  the  weak  convergence  of  the  amplitude  and  the  phase  in  (2.25), 
we  assume  the  following  conditions  hold: 

(2.27)  lim  sup  P{  sup  | Y^(s*)-  Y^(s)|  >  6}  =  0,  V6  >  0 

hTO  e«.[0,e0]  |s'-s|<h  u 

(2.28)  lim  sup  P{  sup  |Y^,(t)  -  Y^(s)]>  6}  =  0,  V6  >  0 

h+0  eftfO.e-J  0<.s,s'<X 

0  | s*-s|<h 

Then  by  invoking  Theorem  1  on  p.  449  and  Theorem  2  on  p.486  in  [4],  we  can 
prove  the  following  lemmas. 

2  £ 

Lemma  2.2  As  e  i  0  and  t  t  «,  with  t  ■  e  t  held  fixed,  we  have  a[Y0(t)]  -=> 
a[CQ(x,t)+  C(i)]  asymptotically  provided  that  the  condition  (2.27))  holds 
and  a  C  C(R)  . 


Lemma  2.3  Suppose  the  condition  (2.28)  holds  and  if  J  g(x)  on  R which  is  posi¬ 
tive,  increasing  as  |x|  -*■  °°  such  that 

sup  sup  E  g[f  (t)]  =  C  <  ® 

E*[0,r-01  tisz.0 


If  b  satisfies 


,  •  lb(t,x) |  _ 

lim  sup  sup  pfxV~  =  0 

t>0  | x | >N  ’ 


then 


<j>(t ,  Ye)  — >  4>(t,Y)  asymptotically  #  . 


Thus  the  above  results  may  be  summarized  as 

Theorem  2.1.  If  the  assumptions  in  Lemmas  2.1  and  2.2  are  satisfied,  then  the 

statistics  of  the  amplitude  a  and  the  phase  $  in  the  solution  (2.25),  in  the 

2 

limit  as  c  1  0,  t  t  *  with  t  *  c  t  fixed,  may  be  computed  according  to  the 


distribution  Y  (t)  given  in  (2.22)  #. 


Remarks : 


(1)  In  general  it  is  not  necessarily  true  that  U =*>  U  =  a(Yp)  exp{ ft , y) } . 

This  is  valid  if  we  can  show  that  a  stronger  convergence,  such  as  the  conver¬ 
gence  in  probability,  takes  place. 

(2)  Theorem  2.1  implies  an  invariance  principle.  For  small  c  and  large  t,  one 
can  approximately  compute  the  statistics  of  a  and  4>  by  means  of  a  Brownian 
motion  (2.22). 

(3)  "a[Yg(t)]  =*>  a[YQ(t)]  asymptotically"  means  that,  writing  Yg(t)  = 

Tg(x,t)  +  £  (t) ,  the  limit  is  taken  with  respect  to  e,  with  £q  and  x  both  held 
fixed.  Similar  interpretation  is  given  for  <t>(t, Y  ). 


3.  RANDOM  HYPERBOLIC  SYSTEMS  IN  ONE  SPACE-DIMENSION 

The  acoustic  wave  propagation  in  one-dimension  may  be  described  by  the 


first-order  system: 


3p  2  3u 

n  3t  Pc  3x  =  ql  ’ 

(o.  1) 

3u  __  3p 

p  3t  3x  q2  ’ 

which  is  subject  to  appropriate  initial  and  boundary  conditions.  In  (3.1), 
p  and  u  denote  the  acoustic  pressure  and  velocity;  p  is  the  mean  density; 
c(t,x,w)  is  the  local  speed  of  sound  which  fluctuates  randomly,  and  q ^ , q2  are 
the  source  terms.  For  convenience,  we  let  =  pu,  =  p  so  that  (3.1)  can 
be  written  as 


3U  .  3U  n 
3t  ’  A  3£  =  Q  • 


where 
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in  the  theory  of  partial  differential  equations  [5],  the  system  (3.2)  is 
usually  analyzed  by  the  method  of  characteristics.  But  this  approach  yields 
a  rather  complicated  statistical  problem.  Such  a  consideration  has  lead  us  to 
a  formal  progressing  wave  solution.  Though  difficult  to  justify  mathematically, 
it  is  physically  appealing  and  has  been  applied  widely  in  studying  time-depen¬ 
dent  wave  propagation  problems.  More  importantly  the  method  may  be  easily 
extended  to  higher  dimensions. 

I.et  us  consider  the  homogeneous  version  of  (3.2): 


(3.4) 


...  3U  .  3U  _  . 

L  U  =  3t  '  A  37=  °  • 


Por  the  first-order  progressing  wave  solution,  we  seek  a  solution  of  the  form 
(p.022,  [51) 


(3.5) 


N 

U(t,x)  -  l  z.  (t ,x)g .  [S(t  ,x) ] 
j=l  3  3 


where  z  is  the  amplitude-vector ;  g  is  the  scalar  wave  form,  and  S  is  the  phase 
function. 

Let 


(3.6) 


U.  =  zjgj(S)  ,  j  =  1,2,... 


Then 


(3.7) 


if  we  choose 


L  ui  =  ®i  Lzi  +  MV'VJ’i 


=  gj  LZ, 


(3.8) 


Mz1  =  (St I  -  Sx A)  zx  =  0 


Let  X  and  z  denote  the  eigenvalue  and  the  right  eigenvector  of  ($XA) , 
respectively.  We  get 


z .  *  p .  z.  , 
1  J  J 


(3.9) 


j  “  1,2,... 
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(3-10)  St  =  A(x,t,Sx)  . 

Thus 

(3.11)  M  =  (U  -  SxA)  , 
and 

(3.12)  L  U2  =  g^Mz2  +  g]Lz2  . 

It  is  easy  to  check  that 

(3-13)  L(U1  +  U2)  =  glLz2 

if  we  choose 

(3.14)  g’  =  gj 

and 


(3.15)  Mz2  =  -Lz1 

This  process  can  be  continued  in  such  a  manner  that 
n 


(3.16) 


U  l  Uj) 

j=l  J 


?lLzn' 


for  any  n  <.  N. 


Wo  tacitly  assumed  that  the  sequence  { zn>  is  decreasing  in  magnitude  so  that 
the  remainder  term  on  the  RHS  of  (3.16)  also  diminishes.  For  simplicity,  we 
stop  at  n  =  2.  In  view  of  (3.9),  (3.15)  becomes 

(3.17)  =  -L(pjZ+) 
which  is  solvable  provided  that 

(3.18)  z"L(p1z+)  =  0  . 

Here  z  denotes  a  left  eigenvector  of  (S^A) .  Hence,  to  a  first-order  approxi¬ 
mation,  the  amptitude  and  the  phase  of  the  wave  function  U  are  determined  by 
(3.10)  and  (3.18). 

According  to  the  Hamilton- Jacobi  theory,  we  introduce  the  "bi-characteristic" 
curve  r  : 

(3.19)  &  =  -  ^  X(s,y,q),  y(t)  =  x, 

3s  *  ly  4W  =  P  =  sx  • 


(3.20) 
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The  curve  I’  =  Y  =  Y(x,s),  0  <.  s  <.  t ,  is  the  physical  characteristics  or  the 
ray.  Note  that  if  X  is  linear  in  q,  then  (3.19)  is  independent  of  (3.20). 
The  amplitude  or  "transport"  equation  (3.18)  can  be  simplified  to  give 

(3.21)  (^  -  X  A)r  +  bp  =  0,  p(0,x)  =  a (x) 

where  p  =  and  X  =  X/S^. 

(3.22)  b  =  (z'l,z+)  with  (z",  z+)  =  1 

Now  the  system  (3. 19) - (3. 21)  resembles  the  system  (2.1)-(2.3)  treated  in 
section  2.  Therefore  the  analysis  there  may  be  carried  over  to  the  present 

case. 


To  illustrate  the  procedure,  we  return  to  the  original  problem  (3.2) 
with  Q  =  0.  It  is  easy  to  see  that  there  are  two  distinct  eigenvalues 

(3.23)  X  =  Xj  2  =  ±  C(t  ,x,w)Sx(t,x,aj) 

The  corresponding  normalized  eigenvector  are 


(3.24) 


'1,2 


_  _1_ 

1 

; 

and  z*  *  -  ! 

1)2  /2 

1 

x:1. 

" 71 

-X,  , 

1,2 

1,2 

The  equation  (3.10),  (3.19)  and  (3.21),  become 

(3.25)  St  =  ±C(t,x,ui)Sx  ,  S(0,x)  =  fc(x)  , 

(3.26)  &  =  ±C(s,y,w),  y  (t)  =  x. 

(3.27)  Pt  ±  Cpx  +  bp  =  0,  p(0,x)  =  a(x) 

where  b  =  (z^  ^z\  2 

Now  suppose  that  the  local  sound  speed  C  fluctuates  about  a  constant  mean 
value  Hq  so  that 

(3.28)  C(t,x,jj)  =  nQ  +  en(t,x,o>) 

where  c  >  0  is  a  small  parameter.  Thus  for  each  characteristic  curve,  we  can 


4  -  “V 
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apply  Theorem  2.1  to  (3.26)  to  obtain  an  asymptotic  limit  law  along  the 

2 

unperturbed  characteristic  curve,  or,  as  f  *  0,  t  t  »  ,  i  »  £  t  fixed, 


(3.29)  Y^  2(s;x,t)  "=s  Y^  2(s;x,t)  asymptotically, 

where  Y(s;x,t)  is  the  Brownian  motion  defined  by  the  Ito  equation  (2.23). 

In  view  of  (3.28),  since  L  is  a  differential  operator,  b  =  0(e).  If 
we  neglect  b  in  (3.27),  the  equations  (3.25)  and  (3.27)  become  identical. 
Their  solutions  are  simply 

(3.30)  S(t,x,u>)  =  0 j  [Yj  (x ,t ,oj) ]  , 

(3.31)  P(t,x,io)  =  a.[Y^(x,t,u))]  ,  Y^(x,t)  =  Y^(0;x,t),  j-1,2  . 

Suppose  we  assume  that  n,  a,  O  satisfy  the  conditions  for  Lemma  2.2.  Then  we 
can  assert  that 

(3.32)  Gj  [Yj  (x,t)]  ->  e^Y  (x,t)] 
and 

(3.33)  a^[Y^(x,t)]  ->  aj[Y^.(x,t)]  ,  j  »  1,2  ,  asymptotically. 

That  is,  for  small  e,  we  may  compute  the  statistics  of  amplitude  and  phase 
for  each  wave  component  by  the  distribution  of  the  Brownian  motion  Y  (0,t). 

To  write  down  pathwise  solution  explicitly,  let  the  initial  state 
11(0, x)  be  written  as 

(3.34)  U(0,x)  *  a1a1(x)z*exp{ik01(x)} 

+a2a2(x)z2exp{ik02(x)}  . 


Then  the  first-order  progressing  wave  approximation  is 

U(t ,x,u>)  ~  ajajfYjix.tilZjexpiikejtYjtx.t)]) 


(3.35) 


+a2a2[Y2(x>t)]z2exp{ik02lY2(x,t)]}  . 


In  view  of  (3.24),  the  eigenvectors  Zj  2  are  random.  But  from  (3.28),  we  have 


(3.36)  z*  =  C  +  0(c),  C  =  — 


+r>  0ex(X±n0t) 


where  C  is  not  random.  Thus  z^  may  be  replaced  by  in  (3.35). 

hot  us  summarize  the  above  results  as  a  theorem. 

Theorem  3.1:  In  the  first-order  progressing  wave  approximation,  the  acoustic 

problem  (3.1)  or  (3.2),  with  q=0  and  the  initial  condition  (3.34)  has  a  sample 

solution  (3.35),  where  may  be  replaced  by  C ^ ,  j=l,2.  Thus,  in  the  diffusion 

2 

limit  v  I  0,  t  t  ®  ,  with  x  =  e  t  fixed,  the  statistics  of  the  amplitude  a. 

and  the  phase  6  ^ ,  j  =  1,2,  can  be  computed  approximately  by  the  probabil  ity 

laws  for  the  Brownian  motions  Y^(x,t)  satisfying 

(3.37)  dY j  =  (-1)-V0dt  +  odw(t) ,  Y^(x,0)  =  x,  j  =  1,2, 

where  w(t)  is  the  standard  Brownian  motion  and  o  is  defined  by  (2.20).  # 


4.  RANDOM  HYPERBOLIC  SYSTEMS  IN  HIGHER  DIMENSIONS 

In  three  dimensions  the  acoustic  system  (2.1)  in  one  dimension  should  be 


modified  to  give 


+  pc"  V*u  =  qj  , 


Here  x  =  (Xj.X2.Xj),  u  *  (Uj.U2.Uj)  and  q2  =  (q2j.q22.923)*  and  they  have 
similar  interpretation  as  before.  For  convenience,  we  set 

(4.2)  >t  .  J  -  1,2,3  . 

Define 


(4.3) 
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(4.4) 


A.  = 
1 


(-1) 


0 

0 

0 

c2s  . 

Ij 


0 

0 

0 


0 

0 

0 


6  =  {*’  1  *  j 
ij  0,  i  t  j,  1,J=1,2, 


(4.5) 


Then  (4.1)  may  be  put  in  the  form 

3 

(4.6)  L  U  =  (3  -  l  A  3.)  U  =  Q  . 

j=l  3  3 


which  is  a  generalization  of  (3,2).  By  an  analogy  to  one-dimensional  case, 
it  is  clear  that  the  progressing  wave  formalism  can  be  extended  to  the  present 
case  in  a  straight-forward  manner.  Therefore  we  shall  only  sketch  the  proce¬ 
dure  without  going  into  the  details. 

In  the  first-order  approximation,  as  before,  let 


(4.7)  U  ~  zg[S(t ,x}]  . 

Then  we  set  Q  =  0  in  (4.6), 

(4.8)  LU  =  gLz  , 
and 

3 

(4.9)  Mz  =  [SJ  -  l  (3  S) A  ]z  -  0  . 

j  =  l  J  3 

The  matrix  M  has  two  repreated  eigenvalues  A  *  0  and  two  distinct  ones 
at 

(4.10)  \1>2  -  iyC  , 


where 


Again  denote  by  z.  (z"  )  the  right  (left)  eigenvector  corresponding  to 

respectively,  with  (z  ,z  )  =1. 

Suppose  that  the  initial  wave  vector  belongs  to  the  subspace  spanned  by 

+  +  . 

Zj,  z^  so  that 

2 

(4.12)  (J(0,x)  =  \  a. a. (x)z+(0,x)exp{ik9 . (x) }  . 

j_1  1  1  J  ) 

Then,  corresponding  to  (3.25)— (3.27) ,  the  amplitude  and  the  phase  are  deter¬ 
mined  by  the  following  set  of  equations 

(4.13)  a  s  =  ±C(t,x,u))|VS|  ,  S(0,x)  =  0(X) 
which  has  the  bi-characteristic  equations: 

(4- 14)  =  _ C(s,y,oi)-j^j— ,  q  t  0,  y(t)  =  x  , 

(4. 15)  =  ±VC(s,y,u)|q|  ,  q(t)  =  p  =  VS  . 

and 

(4. 16)  3tP  T  VS-Vp  +  bp  =  0,  p(0,x)  =  a(x)  . 

In  contrast  with  the  one-dimensional  case,  the  system  of  equations  (4.13)- 

(4.16)  are  highly  nonlinear.  Even  without  randomness,  the  system  is  rather 
difficult  to  deal  with.  In  order  to  proceed,  we  shall  linearize  the  system. 
Again  assume 


(4.17)  C(t,x,u>)  =  n0  +  En(t,x,w). 

Then  (4.15)  implies  q  «  p,  and  (4.14)  yields 


♦ 

s 

r 

k 

(4.18) 

^  =  ±C(s,y,w)p,  y(t)  “  x. 

i 

a 

1 

(4.19) 

P  -  p/Ip!  ■  vs0/ | VSQ |  , 

where  S()(t,x)  denotes  the  phase  function  as  e  =  0.  The  "eikonal"  equation 
(4.13)  may  be  approximated  by 

(4.20)  aS  =  ±  C(t,x,u))p  •  VS,  S(0,x)  =  0 (x)  . 

In  the  same  spirit,  (4.16)  is  replaced  by 

(4.21)  3  p  +  C  p  •  Vp  =  0,  p(0,x)  =  a(x) 

Now  the  linearized  system  (4.18)  -  (4.21)  is  similar  to  (3.25)  -  (3.27)  in 
one-dimension,  and  can  be  treated  accordingly.  We  are  able  to  generalize 
Theorem  3.1  to  the  following. 


Theorem  4.1:  In  the  first-order  progressing  wave  approximation,  the  acoustic 
wave  problem  (4.1)  or  (4.6),  with  Q  =  0  and  the  initial  condition  (4.12),  has 


a  sample  solution  of  the  form 


I  a  a  .[y*(x,t)]exp{ik  0.[y^(x,t)]> 
j-1  3  3  3  3  3 


Further,  in  the  diffusion  limit,  the  statistics  of  the  amplitude  a.  and  the 
phase  0.  can  be  determined  by  the  probability  distribution  of  the  Brownian 
motion  Y^(x,t)  which  satisfies 

(4.23)  dy.  =  ( (-1) juQdt  +  o  dw(t)]pQ  , 


y.(x,0)  =  x,  j  =  1,2, 

where  p_  =  — — ,  w(t)  is  the  standard  BTOwnian  motion,  and  a  is  defined 
°  |V6(x) | 

as  in  (2.20)  #. 


REMARKS 

(1)  The  method  presented  above  can  be  applied  to  other  random  hyperbolic 


systems. 
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(2)  Since,  for  simplicity,  we  assume  the  mean  value  E  C  =  r\  is  constant. 
This  will  limit  the  application  in  underwater  acoustics  to  the  propagation 
over  the  horizontal  range. 

(3)  In  theory  the  straitification  in  depth  may  be  treated  as  well.  The 
difficulty  in  computation  lies  in  the  deterministic,  rather  than  stochastic 
part . 

(4)  As  in  theory  of  geometric  optics,  by  properly  choosing  the  asymptotic 
sequence,  we  should  be  able  to  study  the  statistical  fluctuations  near  a 
smooth  caustics,  which  is  important  in  the  underwater  sound  propagation. 

(5)  Since  the  sample  solutions,  in  the  diffusion  limit,  are  Brownian 
functionals,  the  actual  computation  of  statistical  properties  of  solution 
may  be  carried  out  by  the  method  of  functional  integrals  as  proposed  by  the 
author  [6]. 
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